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We perform a detailed analysis of high-order harmonic generation in diatomic molecules within the 
strong-field approximation, with emphasis on quantum-interference effects. Specifically, we investi- 
gate how the different types of electron orbits, involving one or two centers, affect the interference 
patterns in the spectra. We also briefly address the influence of the choice of gauge, and of the 
initial and final electronic bound states on such patterns. For the length-gauge SFA and undressed 
bound states, there exist additional terms, which can be interpreted as potential energy shifts. If, 
on the one hand, such shifts alter the potential barriers through which the electron initially tun- 
nels, and may lead to a questionable physical interpretation of the features encountered, on the 
other hand they seem to be necessary in order to reproduce the overall maxima and minima in the 
spectra. Indeed, for dressed electronic bound states in the length gauge, or undressed bound states 
in the velocity gauge, for which such shifts are absent, there is a breakdown of the interference 
patterns. In order to avoid such a problem, we provide an alternative pathway for the electron to 
reach the continuum, by means of an additional attosecond-pulse train. A comparison of the purely 
monochromatic case with the situation for which the attosecond pulses are present suggests that 
the patterns are due to the interference between the electron orbits which finish at different centers, 
regardless of whether one or two centers are involved. 
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I. INTRODUCTION 

In the past few years, high-order harmonic genera- 
tion (HHG) and above-threshold ionization (ATI) from 
aligned molecules in strong laser fields of femtosecond 
duration have proven to be a powerful tool for resolving, 
or even controlling, processes in the subfemtosecond and 
subangstrom scale. For instance, one may employ HHG 
and ATI in the tomographic reconstruction of molecular 
orbitals [l], and in the attosecond probing of dynamic 
changes in molecules 0. 

This is possible due to the fact that the physical mecha- 
nisms governing both phenomena take place in a fraction 
of the laser period, i.e., within hundreds of attoseconds 
Q , and involve the recombination or the elastic scatter- 
ing of an electron with its parent molecule [J] . Thereby, 
high-order harmonics or high-energy photoelectrons, re- 
spectively, are generated. Thus, the spectral features are 
highly dependent on the spatial configuration of the ions 
with which the electron rescatters or recombines, and 
yield patterns which are characteristic of the molecule. 
Furthermore, they also depend on the alignment angle of 
the molecule, with respect to the laser-field polarization. 

Due to their simplicity, in particular diatomic 
molecules have been investigated, and minima and max- 
ima have been encountered in their HHG and ATI spec- 
tra. Such patterns have been observed both theoreti- 
cally @, S 0, H, S El HO, El El M El M El E3 

and experimentally [19j, and have been attributed to 
the interference between high- order harmonics or photo- 
electrons generated at different centers in the molecule. 
They are, in a sense, the microscopic counterpart of 
those obtained in a double-slit experiment. Further- 
more, the energy positions of the maxima and min- 



ima depend on the alignment angle and on the inter- 
nuclear distance and, additionally, reflect the bonding 
or antibonding nature of the highest occupied molecu- 
lar orbitals in question [3, El EH . Such features have 
been studied either numerically, by sol ving the time- 
dependent Schrodinger equation (TDSE) pj. El El , or 
semi- analytically, by employing the strong-field approxi- 
mation (SFA) [EEE HI EH El, E El III E3 ■ In Par- 
ticular, in [Tjl, the contributions to the yield from each 
molecular center have been singled out within a TDSE 
computation. Therein, it has been explicitly shown that 
the maxima and minima in the spectra are obtained 
due to the interference between contributions from dif- 
ferent centers, in agreement with the double-slit model 
in [H El. 

Specifically in the SFA framework, the transition am- 
plitude can be written as a multiple integral, with a semi- 
classical action and slowly- varying prefactors. The struc- 
ture of the molecule (and thus its double-slit character) 
can be either incorporated in the prefactors or in the 
action. The latter approach also takes into account pro- 
cesses in which an electron rescatters or recombines with 
a center in the molecule different from the site of its re- 
lease HHEIEIEI. 

An open question is, however, which types of electron 
orbits are responsible for specific interference features. 
For instance, are the dips and maxima originated by the 
interference between orbits in which the electron leaves 
and returns to the same center (regardless of which), or 
between those in which the electron leaves one atom and 
recombines with the other? On the other hand, it could 
also be that the interference patterns result from the com- 
bined effect of all such orbits, and one is not able to at- 
tribute them to specific sets. 



2 



Since the strong-field approximation is a semi- 
analytical method, and allows an immediate association 
with the classical orbits of an electron returning to its 
parent molecule, it appears to be an ideal tool for tack- 
ling this problem. This approximation, however, pos- 
sesses several drawbacks. First, the SFA is gauge de- 
pendent, which leads to different pre-factors and action, 
depending on whether the velocity gauge or the length 
gauge is taken. Second, even for a specific gauge, the pre- 
cise expressions for the pre-factors are not agreed upon, 
and some of them lead to interference patterns which are 
in disagreement with the experiments, and with results 
from other methods [HI, [2(| • 

The main objective of this work is to make a detailed 
assessment of the contributions and relevance of the dif- 
ferent types of recombination scenarios to the above- 
mentioned interference patterns, within the Strong-Field 
Approximation, for a diatomic molecule. Throughout 
the paper, we approximate the strong, infra-red field 
by a linearly polarized monochromatic wave E(t) = 
Eq sin(u!t)e x , and consider a linear combination of atomic 
orbitals (LCAO approximation). In general, we consider 
that the electron reaches the continuum by tunneling 
ionization. As an exception, however, we also take an 
attosecond-pulse train superposed to the monochromatic 
wave [U [H, [HI- The attosecond-pulse train pro- 
vides an additional pathway for the electron to reach the 
continuum, and, recently, has proven to be a powerful 
tool in order to control high-order harmonic generation 
[U, HI and above-threshold ionization [H [13 . In 
the context of the present work, it is a convenient way 
to avoid problems related to spurious potential-energy 
shifts. These shifts are present in the length gauge and 
artificially modify the potential barrier through which the 
electron tunnels. Hence, they may lead to a questionable 
physical interpretation, what the relevance of the differ- 
ent sets of orbits to the patterns concerns. 

This paper is organized as follows. In Sec. HH we 
discuss the strong-field approximation transition ampli- 
tude for high-order harmonic generation, starting from 
the general expressions (Sec. Ill A[) . and, subsequently, 
addressing the specific situation of a diatomic molecule 
in the LCAO approximation (Sec. IIIB[) . Thereby, we 
consider the situation for which the structure of the 
molecule is either incorporated in the pre-factor, or in 
the semiclassical action, in the presence and absence of 
the attosecond-pulse train. When discussing the former 
case, we emphasize the role of the overlap integrals, in 
which the dipole moment and the atomic wave func- 
tions are localized at different centers in the molecule. 
In the latter case, we follow the model in Ref. [ll|, for 
a purely monochromatic driving field, and our previous 
work [231 ] when the attosecond pulses are present, very 
closely. Subsequently (Sec. IIII|) . we investigate the inter- 
ference patterns. Finally, in Sec. IIVI we summarize the 
paper and state our main conclusions. 



II. TRANSITION AMPLITUDES 

A. General expressions 

In general, the strong-field approximation (SFA) con- 
sists in neglecting the influence of the laser field when the 
electron is bound, the atomic or molecular binding poten- 
tial when the electron is in the continuum and the inter- 
nal structure of the system in question, i.e., contributions 
from its excited bound states. The SFA transition am- 
plitude for high-order harmonic generation reads, in the 
specific formulation of Ref. [2(| and in atomic units, 

oo rt 



b Q =1 dt dt' / d 3 Kec(PWMon(p(0) 
J — ocJ — oo J 

exp[iS(t,t', J7,k)] +c.c, 



(1) 



with the action 
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S(t,t',Sl,p)=--J [p + A(T)} 2 dT-I p {t-t') + nt (2) 

and the prefactors d lec (p(t)) = (p(i)| 0dip.e x |^o) an d 
dion(p(0) = (p(i') I Hint{t') \ipo) ■ In the above equations, 
Cdip, e x, H lnt {t'), I p , and fl denote the dipole oper- 
ator, the laser-polarization vector, the interaction with 
the field, the ionization potential, and the harmonic fre- 
quency, respectively. The explicit expressions for p(i) are 
gauge dependent, and will be given below. The above- 
stated equation describes a physical process in which an 
electron, initially in a field- free bound-state {tpo), is cou- 
pled to a Volkov state \p(t')) by the interaction Hi nt (t') 
of the system with the field. Subsequently, it propagates 
in the continuum and is driven back towards its parent 
ion, with which it recombines at a time t, emitting high- 
harmonic radiation of frequency f2. In Eq. {]]), addition- 
ally to the above-mentioned assumptions, the further ap- 
proximation of considering only transitions from a bound 
state to a Volkov state in the dipole moment has been 
made. In the single atom case, this has been justified by 
the fact that the remaining contributions, from the so- 
called "continuum-to-contimmm'' transitions, were very 
small. For a discussion of the various formulations of the 
SFA see, e.g., [H, [H, H3 and in particular [Hj]. 

Due to the fact that it can be carried out almost 
entirely analytically, the SFA is a very powerful ap- 
proach. It possesses, however, the main drawback of 
being gauge-dependent (for general discussions see, e.g., 
Ref. [3(1, and for the specific case of molecules, Ref. 
[ill]). Apart from the obvious fact that the interaction 
Hamiltonians H luii {t'), which are present in di on (p(t')), 
are different in the length and velocity gauges [311 ] . in 
most computations where the SFA is employed, field- 
free bound states are taken, which are not gauge equiva- 
lent. Indeed, a field-free bound state ipg^^ m the length 
gauge would be gauge equivalent to the field-dressed state 
^0 V) ) = Xv*-l V"o L) ) in the velocity gauge, with Xv^i = 



3 



exp[iA(t) • r]. Such a phase shift causes a translation 
p — > p A(t) on a momentum eigenstate |p) . Hence, for 
field-free bound states in both gauges it leads to different 
dipole matrix elements d rcc (p(t)) and di on (p(t)). Explic- 
itly, in the length gauge p(i) = p + A(t), while in the 
velocity gauge p(t) = p. 

For computations involving a single atom, the latter ar- 
tifact can be avoided by placing the system at the origin 
of the coordinate system. For systems composed of sev- 
eral centers, such as molecules, however, this ambiguity 
will always be present. Indeed, in the literature, different 
results have been reported for molecular SFA computa- 
tions in the velocity and in the length gauge @, B El, Q3 • 
In recent, modified versions of the length-gauge SFA, this 
problem has been eliminated for ATI by considering the 



initial bound state 



4 L) ) 



(L) 



Such 



exp[-iA(t')-r] Vo" 
a state is gauge-equivalent to a field-free bound state in 
the velocity gauge, and, physically, may be interpreted 
as a field-dressed state, in which the laser-field polariza- 
tion is included [IT], l32|. For HHG, one may proceed in 
a similar way, with the main difference that the dress- 
ing should also be included in the final state, with which 
the electron recombines. In the dressed modified SFA, 
p(t) = p. 



B. Diatomic molecules 

We will now apply the SFA to a diatomic molecule. 
For this purpose, we will consider the simplest scenario, 
namely a one-electron system and frozen nuclei. Fur- 
thermore, we will assume that the molecular orbital 
from which the electron is released and with which it 
may recombine is a linear combination of atomic or- 
bitals (LCAO approximation). Explicitly, the molecular 
bound-state wave function reads 



Mr) =C^(4 1) (r 1 ) + e4 2) (r 2 )), 



(3) 



where e = ±l,i"i = r — R/2, and r 2 = r + R/2 de- 
note the positions of the centers C\ and C 2 , respectively, 
and Cty is a normalization constant. For homonuclear 

molecules, which we will consider here, tp^ = V'^ = fo- 
The positive and negative signs for e correspond to bond- 
ing and antibonding orbitals, respectively. Within this 
context, there exist two main approaches for comput- 
ing high-order harmonic spectra, which will be discussed 
next. 



1. Pre] "actors 

The simplest and most widely used [f| 0, H, E3) EE 
E3, H3] approach is to incorporate the structure of the 
molecules in the prefactors d ICC (p(t)) and di on (p(*))> 
and to employ the same action S(t,t',Cl, p) as in the 
single-atom case. The multiple integral in the transi- 
tion amplitude JT]) can then either be solved numeri- 



cally, or using saddle point methods [29(. The latter 
procedure can be applied for high enough intensities 
and low enough frequencies, and consists in approximat- 
ing ([T]) by its asymptotic expansion around the coordi- 
nates (t s , t' s , p s ) for which S(t, t' , ft, p) is stationary. This 
implies that d t S{t, t', O, p) = d t ,S(t, t', O, p) = and 
d p S(t, t' , O, p) = 0. In this paper, we employ the specific 
saddle-point approximations discussed in Ref. [331 ]. 

For a single atom placed at the origin of the coordinate 
system, this leads to the equations 



and 



[p + A(t')f = -2I p , 



dr [p + A(t)] = 0, 



2(n - i p ) = [p + A(t)y 



(4) 



(5) 



(6) 



Eq. ((J) expresses the conservation of energy at the time 
t' at which the electron reaches the continuum by tun- 
neling ionization. This equation possesses no real solu- 
tion, which reflects the fact that tunneling has no clas- 
sical counterpart. In the limit I p — > 0, one obtains such 
a condition for a classical particle reaching the contin- 
uum with vanishing drift velocity. Eq. ^ constrains the 
intermediate momentum of the electron, so that it re- 
turns to its parent ion, and, finally, Eq. describes the 
conservation of energy at a later time t, when the elec- 
tron recombines with its parent ion and a high-frequency 
photon of frequency f2 is generated. 

The matrix element d roc (p) = (p ©dip • &x \ipo) then 
reads 



rfrcc(p) 



where 



(2tt) 3 / 2 



e* ,R / 2 I(n) + ee-^- R / 2 I(r 2 )l , (7) 



2"( r i) = J °di P ( ^ | ^ ) ■ e x exp[ip • r^ofa)^, (8) 

and Odi P C 1 2 r2 ) i s the dipole moment. In the length 
gauge, which we are mostly adopting in this paper, 
p(t) = p + A(t). Unless strictly necessary, in order to 
simplify the notation, we do not include the time depen- 
dence in p. The dipole moment can be written in several 
forms. If one considers the length form a natural choice is 
O d jp(r) = — er + er - ! +er 2 = er. Other possibilities are to 
consider the operator Odip in its velocity and acceleration 
forms @, 11. 

Inserting 0di P (r) in Eq. ((8]) yields 

oc +Xj V {vj), with j = 1,2 and v ^ j, (9) 

where 



lj(rj) = J 0dip(>j) • e x cxp[ip • r j ]tp (r j )d 3 r j 



(10) 
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and 



Xj V {rj) = J O dip (r u ) ■ e x exp[ip • r^ofa)^-, v ^ j. 

(11) 

Specifically, if the dipole is in the length form, the above- 
stated integrals read 



(12) 



and 



where 



ZjvOj) = ^H^^ 5 ) + e i^^(P))> (I 3 ) 



(p) 



exp[ip • Tj](po(rj)d 3 rj. 



(14) 



In Xj v , £j = +1 for j = 1 and tj = —1 for j = 2. Eq. 
([7]) is then explicitly written as 



41Up) 



2iCV 
(2tt) 3 / 2 



cos(i?)<9 Px 0(p) 



sin(tf)0(p) 



for bonding molecular orbitals (i.e., e > 0), or 



(15) 



4£(f>) = 



2C</ 



(2tt) 3 / 2 



sin(tf)dpXp) 



cos($)0(p) 



(16) 

in the antibonding case (i.e., e < 0), with $ = p • R/2. 

In Eqs. (fT5|) and (flB|) , the terms with a purely trigono- 
metric dependence on the internuclear distance yield the 
double-slit condition in [18]. The maxima and minima 
in the spectra which are caused by this condition are ex- 
pected to occur for 



p • R = 2nir and p • R = (2n + 1)tt, 



(17) 



respectively, for bonding molecular orbitals (i.e., e > 0). 
For antibonding orbitals, the conditions are reversed, i.e., 
the maxima occur for the odd multiples of 7r and the 
minima for the even multiples. If the velocity gauge is 
taken, the above stated conditions hold for p(t) = p, 
instead of p(t) = p + A(t). This is due to the fact 
that the initial and final free-field sttes are not gauge 
equivalent, as discussed in Sec. IIA. 

The remaining terms grow linearly with the projec- 
tion R x of the internuclear distance along the direction 
of the laser-field polarization, and may lead to unphysi- 
cal results [ll|, [2(| ■ For this reason, they are sometimes 
neglected in the integrals Xj V (rj) There exists, how- 
ever, no rigorous justification for such a procedure. In- 
deed, only recently, it has been shown that such terms 
can be eliminated by considering an additional interac- 
tion which depends on the nuclear coordinate. This in- 
teraction is present in a modified molecular SFA, in its 
dressed and undressed versions fl7| . and leads to contri- 
butions which cancel out the linear term in R x . 



In the length gauge, if the length form of Odip is taken, 
drec(p(t)) = d ion (p(t')), with p(t) = p + A(t), while in 
the velocity gauge, 



^(p)= ^ [ ^f )]2 cos W ^(p), (18) 



or 



4:>(p) = -i ^'pff 12 ^)^), (19) 

with p(t) = p, for bonding and antibonding molecular 
orbitals, respectively. 



2. Modified saddle-point equations 



Physically, if one employs the pre- factors (jl5|) and (|16p , 
this means that one is not modifying the saddle-point 
equations (j4j-([6]). Therefore, the orbits along which the 
electron is moving in the continuum remain the same as 
in the single-center case. This approach is questionable 
in several ways. From the technical viewpoint, there is 
no guarantee that such pre-factors are slowly varying, as 
compared to the semiclassical action, especially for large 
internuclear distances [ill ]. Furthermore, in general, 
they do not incorporate processes in which the electron 
leaves one center of the molecule and recombines with 
the other, which, physically, are expected to be present 
in molecular HHG [llj ]. 

A slightly more sophisticated approach is to exponen- 
tialize the pre-factors obtained in the two-center case and 
incorporate the terms in p • R/2 in the action. This pro- 
cedure has been adopted in [ll| and will be closely fol- 
lowed in this work. For the sake of simplicity, we will 
consider the modified prefactor 



2iC^ 
"(2tt) 3 / 2 



cos I P ' y ) °r, 



(P) 



(20) 



for which the second term in Eq. (fT3|) is absent. Eq. (|20|) 
is also very similar to the dipole matrix element in the 
velocity form, apart from the fact that, in the latter case, 
dpx&ip) is replaced by p^(p) [20(. In the expression for 
the antibonding case, the cosine term in (|20| should be 
replaced by sin(p • R/2). 
This leads to the sum 



(21) 



of the transition amplitudes 

ft 



M 



(2 

x exp[iSj V (p,Q,t,t')}, 



(22) 



with77(p,M') = [dfte<f>(p(t))Y dfa<f>(p(t')). The terms S jv 
correspond to a modified action, which incorporates the 
structure of the molecule. 
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Explicitly, for the undressed length-gauge SFA, 

=S(p,n,t,t') + (-iy+^ 1 (R,t,t') (23) 

and 

s jv = s( P ,n,t,i/) + (-iy+%(R,t,t'), v ?j, (24) 

where £i (R, t,t') = [A(t) - A(t')] ■ R/2 and & [R, t, t') = 
p • R+[A(i) + A(t')} ■ R/2. Eq. ([23]) and ([24]) may be di- 
rectly related to physical processes involving one or two 
centers in the molecule, respectively, as it will be dis- 
cussed next. 

For this purpose, we will solve the multiple integrals 
in ([22]) employing saddle-point methods. The condi- 
tions d p S jv (p,Sl,t,t') = 0, d t S jv (p,Sl,t,t) = and 
dt>Sjv(p,tl,t,t') — upon the derivative of the action 
yield saddle-point equations, which, as in the previous 
section, can be related to the orbits of an electron re- 
combining with the molecule. Explicitly, for the modified 
action Sjj , [Eq. ([23))] , the saddle-point equations read 



A(t')] 2 



for j = 1, or 



[p + Mt')} 2 



■A(t)f 



= -I p - E(0 • R/2, (25) 



Sl-I p + E(i) • R/2, (26) 



-I p + E(i') ■ R/2, (27) 



9.-I p - E(i) • R/2, (28) 



for j = 2. The saddle point equations ([23)) and ([2"7| 
correspond to the tunnel ionization process from cen- 
ter C\ and C2, respectively. Curiously, both equations 
contain extra terms, if compared to equation ^ for the 
single-atom case. Such terms are dependent on the inter- 
nuclear distance R and the external laser field E(i') at 
the time the electron is freed, and may be interpreted as 
potential-energy shifts in the barrier through which the 
electron tunnels out. Similar terms are also observed in 
Eqs. ([2"6"]) and (|28[) for the energy conservation at the time 
the electron recombines, as compared to the single-atom 
expression ©. 

The remaining saddle point equation is given by the 
same expression as in the single-atom case, i.e., Eq. ([3]), 
and means that, for M\\ and M22, the electron is ejected 
and returns to the same center. If on the other hand, we 
consider the modified action ([2~4"]) , this yields 



(29) 



/ [p + A(s)]ds ± R = 0, 
Jv 



and positive signs refer to the transition amplitudes M12 
(center C\ to center C2) and M 2 \ (center C2 to center 
Ci), respectively. In the former case, the remaining sad- 
dle points are given by ([2"5| and (|28[) . i.e., the electron 
tunnels from C\ and recombines with C'2 , whereas in the 
latter case they are given by ([2"T]) and ([2"6")1 , which physi- 
cally, expresses the fact that the electron is ejected at C2 
and recombines with C\. 

The energy shifts ±E(r) - R/2, T = t,t' in ([25 ]) -([28 )1 are 
absent in the velocity gauge SFA [Tlj, and in a modified 
length-gauge SFA, in which the electric field polarization 
is incorporated in the initial and final electronic bound 
states. In both cases, the action Sjj, associated to orbits 
involving only one center, is given by the single-atom 
expression |2]), which leads to the saddle-point equations 
(ID)-©. The action Sj V related to two-center orbits reads 



s jv = s(p, n, t, f) + (-ir +1 P ■ r, v ± j 



(30) 



which, physically, mean that the electron is leaving from 
one center and recombining with the other. The negative 



The above-stated expression leads to the single- 
atom equations ([!]), ([6]) for tunneling and rescatter- 
ing, together with the two-center return condition 
([29)) . The prefactors, however, are different in both 
cases. In the dressed modified length-gauge SFA, 
i] D (p,t,t') = [dp 4>(p)]* dp 4>(p), while in the velocity 
gauge V V (p,t,t') = [9 p ^(p)]^(p)[p + A(t')] 2 /2- 

One should note that, within the specific model em- 
ployed here, the physical process in which the electron 
moves directly from one center to the other, without 
reaching the continuum, is not being considered. Such 
a process leads to a strong set of harmonics in the low- 
energy range of the spectra. Since, however, we are fo- 
cusing on the plateau harmonics, the contributions from 
this extremely short set of orbits are not of interest to 
the present discussion. For a detailed study of this case, 
see, e.g., 



3. Additional attosecond pulses 

The role of the energy shifts observed in ([25 )) -([28 )) is not 
well understood. A way of eliminating such terms is to 
modify the length-gauge SFA and include the influence 
of the laser field in the initial and final states. How- 
ever, even without such modifications, it is possible to 
provide an additional pathway for the electron to reach 
the continuum, so that, at least in the context of tunnel- 
ing ionization, these shifts can be avoided. For instance, 
if the electron is ejected by a high-frequency photon, it 
does not have to tunnel through potential barriers with 
energy shifts whose physical meaning is not clear. Such a 
pathway can be provided by a time-delayed attosecond- 
pulse train Eh(f) superposed to a strong, near infra-red 
field E; (t) = Eq sin u>te x . Indeed, it has been recently 
shown that such pulses can be used to control the elec- 
tron ejection in the continuum, and thus high-harmonic 
generation and above-threshold ionization [2]], • 

In [H, |24j, we employed such a scheme to control 
quantum-interference effects for high-harmonic genera- 
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tion and above-threshold ionization for the single-atom 
case, within the SFA framework. Our previous findings 
suggest that the probability of the electron reaching the 
continuum, in case it is ejected by the attosecond pulses, 
is roughly the same for all sets of orbits. Indeed, it ap- 
pears that the sole, or at least main factor determining 
the intensities in the spectra is the excursion time of the 
electron in the continuum. In the specific case studied 
in [23j|, there was a set of very short orbits, which led to 
particularly strong harmonics. Therefore, an attosecond 
pulse train superposed to a strong laser field is an ideal 
setup to avoid any artifacts due to modified tunneling 
conditions. 

In [23|], we have approximated the attosecond pulse 
train by a sum of Dirac-Delta functions in the time do- 
main. This yields 



In case one considers the transition amplitudes (|22| , this 
yields 



M ju = 



■ f t 771 00 /* + 00 r- 

|^E(- 1 )7_ o f/^cp,*,0(33) 



x exp[iSj U (p,Q, t,t') 



(34) 



In both equations, t' = t e i+mt /to. The saddle-point equa- 
tions ([5]) and ([6]) for the single atom case can then be 
combined as 



sin ojt — (—1)™ sintutd 



[ui(t - t d ) - rm] I cos ut t J ^ v 



(35) 



Eh(t) = E h ir 



(~l) n r . rar 

y —J-5(t t d )e a 

a(t) to 



(31) 



where to, Eh, and <r(t) denote the laser field frequency, 
the attosecond-pulse strength and the train temporal en- 
velope, respectively. This approximation is the limiting 
case for a train composed of an infinite set of harmon- 
ics, and has the main advantage of allowing an analytic 
treatment of the transition amplitudes involved up to at 
most one numerical integration. Furthermore, it is a rea- 
sonable asymptotic limit for pulses composed by a large 
high-harmonic set. We consider here a = const, which, 
physically, corresponds to an infinitely long attosecond- 
pulse train. Clearly, the total field is given by E(i) = 
E ; (f)+E h (t). 

We will now assume that the attosecond pulses are the 
only cause of ionization and that the subsequent propaga- 
tion of the electron in the continuum is determined only 
by the monochromatic field. Hence, E(i') ~ Eh(t') and 
A(t) ~ Ai(t) in Eqs. (T]) and (J22J)- This eliminates the 
integral over the ionization time in the transition ampli- 
tudes. Hence, the values for which the single atom-action 
S(t, t', SI, k), or the modified action Sj U is stationary must 
be determined only with respect to the variables t and 
p. Physically this means that the recombination and 
return conditions remain the same, with regard to the 
purely monochromatic case, but that there is no longer 
a saddle-point equation constraining the initial electron 
momentum. In fact, the electron is being ejected in the 
continuum with any of the energies Nto — I p , since all 
harmonics composing the train are equivalent. For the 
other extreme limit, namely a high-frequency monochro- 



matic wave, we refer to 241 ]. where we provide a detailed 



discussion within the SFA. 

Explicitly, if the action is not modified, these assump- 
tions lead to the transition amplitude 



which will give the return times t. If, on the other hand, 
these equations are modified, it is possible to distinguish 
four main scenarios. Specifically, for the processes in 
which the electron leaves and returns to the same center, 
the saddle-point equations differ from (f3"5)) only in a shift 



Ip — * Ip ± E/(<) • R/2. The negative and positive signs 
correspond to Mxx and M22, respectively. For the sce- 
narios involving two centers, the saddle-point equations 
read 



smujt — (—1)" smustd + e 



2U V 



[uj(t — td) — mr] cos tut =F 



I SI 



2U V 



(36) 



with I p = I p 



eEi(t) ■ R/2. The case e = -1 and e = 
corresponds to Mxi and M21 , respectively. 



III. HARMONIC SPECTRA 

We will now present high-harmonic spectra, in the 
presence and absence of the attosecond pulses. We re- 
strict the electron ionization times to the first half cycle 
of the driving field. We also consider the six shortest pairs 
of orbits for the returning electron. Due to wave-packet 
spreading, the contributions from the remaining pairs are 
negligible [35| . For simplicity, we employ a bonding com- 
bination of Is states, for which 



(P) 



1 



[P 2 



(37) 



(D) _ inC^Eh ^ 

CT ( 27r )3/2 

><Cc(pWMon(p(0) 



and assume that the molecule is aligned parallel to the 
laser-field polarization, so that R x = R. 



Prefactors 



dt / d 3 P exp[iS(p,n,t,t')} 

v ' We will commence by considering the single-center ac- 

(32) tion @ and the prefactors discussed in the previous sec- 
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FIG. 1: (Color Online) High-harmonic spectra computed em- 
ploying the single-atom orbits and two center prefactors, us- 
ing a bonding (e = +1), linear combination of atomic orbitals, 
for internuclear distances 1 a.u.< R < 8 a.u. The upper and 
lower panels correspond to the dipole prefactors (I15p and (I20|) . 
respectively. The atomic system was chosen as H^, which was 
approximated by the linear combination of Is atomic orbitals 
with I p = 0.5 a.u.. We took the the driving field intensity 
and frequency to be I — 1 x 10 15 W/cm 2 , and uj = 0.057 a.u., 
respectively. 

tion. In Fig. [T] we depict high-harmonic spectra com- 
puted for a wide range of internuclear distances, employ- 
ing the prefactor ([Pa]) or the modified expression PD|) , for 
which the linear term in R x is absent [upper and lower 
panel, respectively]. The figure illustrates how the linear 
term masks the interference patterns. In fact, for a.u. 
1 < R < 3 a.u., it counterbalances the influence of the 
purely trigonometric term, and no clear minima and max- 
ima are observed. As the internuclear distance increases, 
this term starts to play a dominant role, displacing the 
maxima and minima away from the double-slit condition 
(I17p . Furthermore, it seems that the patterns become 
more distinct with increasing harmonic order. 

A rough estimate of the influence of each term in (fT5|) 
on the spectra agrees with Fig. 1. Since the trigono- 
metric functions in (|15[) are bounded, the ratio be- 
tween the maxima caused by each term will be ? = 
R\ cos0</>(p)/(29 Pa .0(p))|, where 6 is the alignment angle 
and </>(p) is given by Eq. (|14p. If <; ~ 1, the maxima will 
possess the same order of magnitude and there will be 
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FIG. 2: (Color online) Spectra computed employing the 
single-atom orbits and two center prefactors (black thin lines), 
as compared to those obtained employing modified saddle- 
point equations (green thick lines). We consider here the 
modified length form (|20|l of the dipole operator, which ex- 
cludes the term with a linear dependence on R x . The atomic 
system was chosen as which was approximated by the 
linear combination of Is atomic orbitals with I p = 0.5 a.u.. 
The internuclear distance and the alignment angle are R — 5 
a.u., and 8 = 0, respectively. The driving field intensity and 
frequency are given by I — 1 x 10 15 W/cm 2 , and uj = 0.057 
a.u., respectively. The interference minima are indicated by 
the arrows in the figure. 

no noticeable modulation, while if q < 1 the double-slit 
physical picture may still be reproduced. For q > 1, how- 
ever, one expects that the linear term in R will prevail. 
Hence, the critical value for the internuclear distance is 
R c = 2| sec#d Px 0(p)/</>(p)|. This expression depends on 
the bound states with which the electron recombines, and 
also on the harmonic energy. For instance, specifically for 
Is states, R c ~ 4p x sec#/(p 2 + 2I p ). Above the ioniza- 
tion threshold, according to Eq. ©, R c ~ 2p x sec9/rt. 
Hence, one expects the linear term to be more promi- 
nent as the harmonic energy increases, leading to clearer, 
though incorrect, patterns. In order to avoid such prob- 
lems, we will employ the prefactor (|20|) throughout. 

In Fig. [2l we compare spectra computed using either 
the pre- factor (fSTTf) or modified saddle-point equations. 
Both spectra are very similar, with maxima and min- 
ima at harmonic frequencies O = I p + n 2 n 2 /(2R 2 ), as 
expected from the double-slit condition. This similarity 
holds not only for the gross features, but, additionally, for 
the substructure caused by other types of quantum inter- 
ference. Close to the minima, however, the yield from the 
latter case is larger. Nevertheless, the very good overall 
agreement shows that, in fact, the patterns obtained can 
also be interpreted as the result of the quantum interfer- 
ence between different types of electron orbits. 
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B. Interference effects 

In order, however, to identify which sets of orbits cause 
the dips and the maxima, we will analyze the interference 
between their individual contributions. Such results are 
shown in Fig. [3] In the upper panels, we present the 
spectra computed from topologically similar sets of or- 
bits, i.e., from processes involving only one, or two cen- 
ters [panel [3] (a)]. In this case, the main interference pat- 
terns are absent. This strongly suggests that they are 
due to the quantum interference of topologically different 
sets of orbits: the orbits along which an electron leaves 
and returns to the same center, and those along which it 
reaches the continuum at one center and recombines with 
the other. Physically, this could be attributed to the fact 
that, in this case, there would be an appreciable phase 
difference between the two sets of orbits, since the latter 
are much longer than the former. This phase difference 
would cause the overall modulation. 

Hence, there are two remaining possibilities. Con- 
cretely, the modulation can be due to the quantum inter- 
ference either between processes in which the electron 
leaves from different centers and recombines with the 
same center (i.e., between the orbits which start at Cj 
and end at C„, with v ^ j, and those starting and end- 
ing at C„, and j = 1,2), or between those in which an 
electron starts at the same center and recombines with 
different centers (i.e., the electron is ejected at Cj and 
recombines at C„, v ^ j, or it is freed and recombines at 
Cj, with j = 1, 2). In panel [21(b), we consider the former 
processes, whereas in panel [3] (c) we depict the latter. 
Interestingly, only in case the electron leaves from the 
same center, the interference patterns are present. Fur- 
thermore, there is a difference in roughly four orders of 
magnitude between the two types of contributions. Such 
a difference is absent in the other cases. Additionally, 
the full contributions to the yield are practically indistin- 
guishable from the transition probability \AI21 + M22I 2 ! , 
from the orbits starting at C%. 

In order to understand this better, one must have a 
closer look at the individual contributions from different 
sets of orbits, and, in particular their orders of magni- 
tude. Such contributions, depicted in Panels (a) and (b) 
of Fig. 21 show that the transition probabilities I-M22I 2 
and |M2i| 2 , which correspond to the orbits starting from 
the center C2 , are roughly four orders of magnitude larger 
than I-M12I 2 and |Mn| 2 , i.e., than those from the orbits 
starting at C\. Therefore, it is not surprising that the 
yield is dominated by \M%\ + M22I 2 in the previous fig- 
ure. Furthermore, these results exhibit no maxima and 
minima. Hence, they support the assumption that such 
features are due to the interference of different types of 
orbits. Finally, the contributions from orbits starting 
at the same center possess the same order of magnitude. 
This suggests that tunneling ionization is the main mech- 
anism determining the relevance of a particular type of 
orbits to the spectra, and that are local differences in the 
barrier through which the electron must tunnel, depend- 
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FIG. 3: (Color Online) Contributions to the high-harmonic 
yield from the quantum interference between different types 
of orbits, for internuclear distance R — 5 a.u. The remaining 
parameters are the same as in Fig. 3. Panel (a): Orbits 
involving similar scattering scenarios, i. e., \Mu +M22I 2 , and 
\Mi2 + M21I 2 . Panel (b): Orbits ending at the same center, 
i.e., |Mii+M 2 i| 2 and A/12 + M22I 2 . Panel (e): Orbits starting 
at the same center, i.e., [Mn + M 12 \ 2 and \M 2 i + M 22 | 2 . For 
comparison, the full contributions \Mn + M22 + Mil + M12I 2 
are displayed as the light gray circles in the picture. 



ing on the center it starts from. 

An inspection of the imaginary parts Im[i'] of the start 
times provides additional insight into this problem. Due 
to the fact that tunneling ionization is a process which 
has no classical counterpart, this quantity is always non- 
vanishing, even if the energy range in question is lower 
than the maximal harmonic energy. The larger Im[i'[ is, 
the less probable it will be that tunneling ionization takes 
place. Such an interpretation has been successfully em- 
ployed in (36j in order to determine the dominant pairs of 
orbits, in the context of nonsequential double ionization 
with few-cycle laser pulses, and will be also considered 
in this work. For that purpose, we will take the shortest 
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FIG. 4: (Color online) Contributions to the high-harmonic 
yield from specific types of orbits, for internuclear distance 
R — 5 a.u. The remaining parameters are the same as in 
Fig. 3. Panel (a): Yield from the orbits starting and ending 
at the same center, i. e., transition probabilities |A4n| 2 and 
IM22I 2 . Panel (b): Yield from the orbits starting and ending 
at different centers, i.e., transition probabilities |Mi2| 2 and 
|M 21 | 2 . 



pairs of orbits utilized in the computation of the tran- 
sition probabilities in Fig. [4] and, for each case, display 
Im[i'] . These are the dominant pairs of orbits. The longer 
pairs have a less significant influence on the spectra, due 
to the spreading of the electronic wave packet [35| . 

Such results are depicted in Fig. [5] Clearly, Im[i'] 
is around four times larger for the orbits starting from 
the center C\, as compared to those starting from GV 
This means that, in order to reach the continuum, the 
electron must overcome a larger barrier if it comes from 
C\. Since, roughly speaking, the ionization probability 
per unit time decreases exponentially with Im[£'], one 
expects the contributions from the orbits starting from 
C\ to be around four orders of magnitude smaller than 
those from the orbits starting at C2. An inspection of 
Fig. [4] shows that this is indeed the case. 

The above-stated effect could, however, be an artifact 
of the strong-field approximation in the length gauge. 
Indeed, the terms ±E(i') • R in the saddle-point equa- 
tions (|25p and (P?)) can be interpreted as potential en- 
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FIG. 5: (Color online) Imaginary parts of the start times t' 
of the shortest pairs of orbits, which contribute to the matrix 
elements M22, and M12, [panel (a)], and to M21, and Mn 
[panel (b)], for the same field and atomic parameters as in 
the previous figure. 



ergy shifts, due to the fact that the electron is displaced 
from the origin [Tl|, E3 • Such terms increase or sink the 
potential barrier for C\ or C2, respectively, and, conse- 
quently, change the orders of magnitude in the contribu- 
tions starting from different centers. Even though, as a 
whole, the results match those obtained by other means, 
their physical interpretation is controversial. 

One may, however, avoid this problem by providing 
an additional pathway for the electron to reach the con- 
tinuum. For that purpose, we shall superpose a time- 
delayed attosecond pulse train to the strong laser field, 
employing the model discussed in Sec. Ill B 31 and in our 
previous work [23L l24j . The maximal harmonic ener- 
gies for this specific model are strongly dependent on 
the time delay td between the attosecond-pulse train the 
infra-red field, extending from the ionization potential, 
for (f> = 0.75-7T to I p + 1.8Up, for <fi = mr. Furthermore, 
there exist many intermediate delays, for which a double 
plateau is present. This substructure may be detrimen- 
tal to the identification and physical interpretation of the 
interference patterns. In order to avoid such problems, 
we will consider here vanishing time delay, i.e., </> = 0. 

In Fig. we depict the spectra obtained for a diatomic 
molecule subjected to such a field, assuming either a two- 
center prefactor and the single-center saddle-point equa- 
tion (|35p . or the modified saddle point equations 
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FIG. 6: (Color Online) High-harmonic spectra for a molecule 
aligned parallel to the laser-field polarization, and internu- 
clear distance R = 3.5 a.u, for an attosecond pulse train 
superposed to a strong, near-infrared field, and exhibiting 
vanishing time-delay (cj> = 0) with respect to it. We con- 
sider a bonding combination of Is atomic orbitals, and take 
the intensity of the attosecond-pulse train to be /j, = J;/ 10. 
The remaining parameters are the same as in the previous 
figures. Panel (a): spectra computed employing the modi- 
fied version of the two-center prefactors in the length form 
[Eq. (|2Up ] and single-atom saddle-point equations, compared 
to that obtained using modified saddle-point equations. Panel 
(b): Contributions from individual scattering scenarios, i.e., 
from |A/jv| 2 , with j = 1,2 and v — 1,2. The interference 
minima are indicated by arrows in the figure. 



[Fig. [HI (a)]. Both computations exhibit a minimum near 
the harmonic frequency £1 — 71a;, in agreement with Eq. 
(fTTj) . If only the contributions \Mj V \ 2 from the individual 
scattering scenarios are taken, such a minimum is absent 
[Fig. [H(b)]. Therefore, it is due to interference effects 
between different sets of orbits. One should note that, in 
contrast to the purely monochromatic case, all contribu- 
tions exhibit the same order of magnitude. This is due 
to the fact that, if the attosecond pulses are present, the 
electron is being ejected in the continuum with roughly 
the same probability, regardless of the center it left from. 

The precise role of the various recombination scenar- 
ios is illustrated in Fig. [7) For clarity, we concentrate 
on the plateau region around the interference minimum. 
The main difference observed, with regard to the purely 
monochromatic case, is that the overall shape of the spec- 
trum, and consequently its minimum, is due to the collec- 



FIG. 7: (Color Online) Contributions to the high-harmonic 
yield from the quantum interference between different types 
of orbits, for the same field and molecular parameters as in 
Fig. 7. Panel (a): Contributions from topologically similar 
scattering scenarios, i. e., contributions from [Mil + M22I 2 
and I-M12 + M21I 2 . Panel (b): Contributions from orbits start- 
ing at the same center. Panel (c): Contributions from orbits 
ending at the same center. For comparison, the overall spectra 
are displayed as the blue symbols in the figure. The interfer- 
ence minimum is indicated by the vertical lines near fi = 71a;. 



tive interference of several types of orbits. This is in con- 
trast to the previous results, for which they were caused 
by the processes starting at a center Cj and ending at 
different centers, i.e., \Mj v + Mjj\ 2 , with v 7^ j and v,j 
= (1,2). Indeed, a modulation is even present for the 
probability \M\\ + M22I 2 involving only one-center sce- 
narios. This is shown in Fig. [JJ (a), and contradicts the 
previously made assumption that such features are due 
to the interference between topologically different sets of 
orbits. In fact, it seems that the absence of overall max- 
ima and minima for the one-center contributions, in the 
purely monochromatic case [Fig. [31(a)], is due to the dif- 
ferent orders of magnitude for Mix an d M22 [c.f. Fig. 
11(a)]. 

Furthermore, one needs several different processes in 
order to obtain the correct position of the minimum. For 
instance, in Fig. [71(a), the overall spectrum closely fol- 
lows [Mil +M22I 2 in the low-energy region. In the vicin- 
ity of the minimum and for higher energies, however, it 
follows neither such contrbutions nor \M\z 4- Af2i| 2 , from 
two-center processes. In Fig. [71(b), where we display the 
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contributions \Mjj + Mj U \ 2 , with v =^ j and v,j — (1, 2), 
from the orbits starting from the same center, the spec- 
trum closely follows I-M21 + M22I 2 before the minimum, 
and \M\i + Mi2\ 2 after the minimum. 

The remaining panel [Fig. [3(c)] depicts the contribu- 
tions from \M UV + Mj V \ 2 , with v 7^ j and u, j = (1,2), 
which give the orbits finishing at the same center. In this 
case, the interference minimum is absent. This is a strong 
evidence that the relevant condition for the presence of 
such features is that the orbits taken into account end at 
different centers, instead of being topologically different 
(which is the case for both Fig. [3(b) and Fig. 0(c)). 
Therefore, these results agree with the double-slit pic- 
ture, which has been put across in (l8l |. 

The additional attosecond pulses have the advantage 
of not introducing changes in the standard length-gauge 
SFA formulation. They modify, however, the physics of 
the problem, since they provide a different mechanism 
for the electron to reach the continuum. Clearly, there is 
also the possibility of eliminating the spurious potential 
energy shifts, by considering a different version of the 
SFA, such as the field-dressed length-gauge formulation 
proposed in [13], or the velocity- gauge formulation. 

In Fig. [51 we display high-order harmonic spectra for 
the field-dressed SFA in the length gauge, and for the 
field-undressed SFA in the velocity gauge [panels (a) and 
(b), respectively]. For simplicity, we exhibit the results 
obtained with the modified prefactor ([20)l . instead of a 
modified action. We also provide curves for which only 
the cosine term has been set to one, in order to facilitate 
the identification of interference effects. As an overall 
feature, we do not observe the interference patterns ex- 
hibited in the previous figures. Indeed, the curves with 
and without the cosine term are very similar. There is, 
however, a minimum near f2 = 1150-1 for the former case. 

For the dressed length-gauge SFA, the above-stated 
features can be attributed to the fact that the condition 
for maxima and minima is now given by (|17|) . with p 
instead of p = p + A(t). The harmonic frequencies for 
which they occur can be easily obtained from condition 
([6|), and are given by 

Q = I p + [n 2 n 2 /Rl + 2rmA{t) / R x + A 2 (*)] /2. (38) 

An upper bound for fl can be estimated as follows. At 
the electron return times, the vector potential is roughly 
^CO ~ 2-\/Up- This yields, for the parameters in Fig. [SJ 
ft ~ 121a;, which is slightly larger than the minimum 
encountered. 

A breakdown of the interference patterns also occurs 
in the velocity gauge, for the very same reasons. In- 
deed, the interference condition for the SFA in the ve- 
locity gauge and for the field-dressed SFA in the length 
gauge are identical. This is a direct consequence of 
the fact that field-free initial and final electron states 
in the velocity gauge are gauge-equivalent to the field- 
dressed states considered in this paper. This gauge equiv- 
alence will lead to identical recombination form factors 
^roc(p) = (p| Cdip-Gx \ipo) . Since the interference condi- 
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FIG. 8: (Color Online) High-order harmonic spectra for the 
same field and molecular parameters as in Fig. 2, but com- 
puted with the field-dressed modified length-gauge formula- 
tion of the strong-field approximation [Panel (a)], compared 
to its field- undressed velocity-gauge counterpart [Panel (b)]. 
The arrows in the figure indicate the harmonic order for which 
a minimum is observed, and the vertical line marks the rough 
estimate for such a minimum. The solid and dotted lines cor- 
respond to the prefactor (19) and to the situation for which 
cosp ■ R/2 has been set to one, respectively. 



tions (|17p are mainly determined by d rcc (p), they will 
be the same. Hence, the harmonic orders for which 
the maxima and minima occur are given by Eq. I|38[). 
and therefore are unrealistically high. The discrepan- 
cies between both yields stem from the form factors 
^ion(p) = (p| Hi n t(t') \ipo), which are gauge-dependent. 

One may, however, consider field-dressed states in the 
velocity gauge, which are gauge-equivalent to field-free 
states in the length gauge. This is achieved by apply- 
ing the transformation Xv*-i — exp[iA(i) -r] in the initial 
and final length-gauge electronic bound states. This leads 
to a shift p — > p + A(i) on a momentum eigenstatc |p), 
which is exactly the opposite shift induced in the field- 
dressed length-gauge SFA. This shift has the main con- 
sequence that the interference condition (|17p now holds 
for p = p + A(t), even in the velocity gauge. The results 
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FIG. 9: (Color Online) High-order harmonic spectra for the 
same field and molecular parameters as in Fig. 2, but com- 
puted with the field-dressed modified velocity-gauge formula- 
tion of the strong-field approximation. In order to facilitate 
the comparison, the undressed length-gauge version of the 
SFA is provided, and the velocity-gauge yield has been nor- 
malized in approximately 2 orders of magnitude, to match the 
length-gauge results. In both cases, for simplicity, the mod- 
ified prefactor (18) has been employed, instead of modified 
saddle-point equations. 



obtained employing such dressed states in the velocity- 
gauge SFA, depicted in Fig. 9, are indeed very similar to 
those obtained using its field undressed, length gauge for- 
mulation. In fact, we have observed mainly quantitative 
differences, due to different prefactors c?i on (p) [s3|- 

One should note, however, that the transformation 
Xv^i introduces the same additional potential energy 
shifts in the modified action as in the undressed length- 
gauge case. Thus, the price one pays for recovering the 
correct interference conditions is the loss of a direct con- 
nection to simple classical models. 



IV. CONCLUSIONS 

The present results support the viewpoint that the 
maxima and minima in the high-order harmonic spectra 
of diatomic molecules are due to the interference of elec- 
tron orbits finishing at different centers in the molecule. 
This seems to hold regardless of whether the electron has 
been released in one center at the molecule and recom- 
bines with a different center, or whether it is ejected and 
returns to the same center. Such conclusions have been 
reached by employing modified saddle-point equations, 
within the strong-field approximation. These modifica- 
tions lead to orbits involving different centers, and are a 
slightly more refined approach than the standard proce- 
dure of considering single-center saddle-point equations 



and modified prefactors. 

In this framework, we considered that the electron has 
been ejected by tunneling ionization and by an additional 
attosecond-pulse train, and compared the similarities and 
differences from both physical situations. In the former 
case, depending from which center the electron is leav- 
ing, it must overcome unequal potential barriers to reach 
the continuum, whereas in the latter case it is ejected 
with roughly equal probability. Especially in the purely 
monochromatic case, we observed that the contributions 
from orbits starting at one of the centers, namely C2, 
were much larger than those from C\ , due to a narrower 
potential barrier. In particular, there exist potential- 
energy shifts which are proportional to the electric field at 
the ionization time and the internuclear distance, which 
cause such differences. It is however noteworthy that the 
electron excursion times have been confined to the first 
half-cycle of the laser field. If the other half-cycle had 
been taken, the potential barrier would reverse and the 
contributions starting from C\ would be more prominent. 

The above-mentioned potential-energy shifts, however, 
do not possess a clear-cut physical interpretation. In fact, 
they are only present in the standard length-gauge formu- 
lation of the strong-field approximation, i.e., if field- free 
bound states are taken when the electron is ejected and 
recombines, and are the source of several problems. For 
instance, they reflect the fact that the SFA is translation- 
dependent. Moreover, due to their existence, it is diffi- 
cult to establish an immediate connection between this 
approach and the classical equations of motion of an elec- 
tron in the laser field. 

On the other hand, it seems that such shifts are nec- 
essary in order to obtain the correct energy position of 
the maxima and minima in the HHG spectra. In fact, an 
improved formulation of the SFA, in which the influence 
of the laser field is included in the electron bound states, 
restores its translation invariance, provides an unprob- 
lematic classical limit (l7j . but yields incorrect energy 
positions for the interference patterns. This discrepancy 
is related to the fact that the field dressing alters the 
double-slit interference conditions (fT7|) . A similar ab- 
sence of interference features has been reported very re- 
cently in Ref. (38|, for HHG computations using a field- 
dressed version of the SFA in the length gauge. 

Finally, when employing different gauges, from our re- 
sults it is clear that the dressing of the initial and the 
final states plays a far more important role than the dif- 
ferent interaction Hamiltonians Hi nt (t'). If the dressing 
is applied consistently so that the the electronic bound 
states are gauge-equivalent, the interference patterns will 
remain the same. This is due to the fact that the inter- 
ference condition (|17| will then remain invariant. For 
instance, the spectra obtained in the velocity- gauge SFA 
with undressed states are very similar to those computed 
in the length gauge with field dressed states. The same 
holds for the undressed length-gauge, and the dressed 
velocity-gauge SFA spectra. In the two former cases, 
there is a breakdown of the interference patterns, as com- 
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pared to the field-undressed length gauge SFA. However, 
such patterns can be restored in the velocity gauge, by 
dressing the electronic bound states appropriately. 
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